Abstract. We address the question of the dynamics in 
the inner 50 pc of the Galactic Centre. In a first step we 
investigate the cloud-cloud collision rate in the Circum- 
nuclear Disk (CND) with the help of a three dimensional 
N-body code using gas particles that can have inelastic 
collisions. The CND might be a longer lived structure than 
previously assumed. The whole disk-like structure of the 
CND can thus survive for several million years. A realis- 
tic simulation of the CND shows the observed disk height 
' ' structure. In a second step the environment of the CND is 
. taken into account. Retrograde and prograde encounters of 
' a cloud of several 10^ Mq falling onto an already existing 
, nuclear disk using different energy loss rates per collision 
5h ' are simulated. The influence of the energy loss rate per 
collision on the evolution of the mass accretion and cloud 
collision rates is strongest for a prograde encounter. A 
(y-j . composite data cube of two different snapshots of a pro- 
^SJ ' grade encounter together with the CND shows striking 
similarity with the observed Sgr A cloud complex. The 
current appearance of the Galactic Centre environment 
can thus be explained by at least two dynamically dis- 
' tinct features together with the CND. The current mass 
■r^lj- ,: accretion rate within the CND ranges between 10~^ and 
lO"** M0yr~^. It can rise up to several 10^^ M0yr^^ 
during massive accretion events. 
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1. Introduction 

The Galactic Centre region is an ideal laboratory for 
studying the fueling of a central black hole in great detail. 
Eckart & Genzel (1996) have observed the proper motions 
of bright stars in the central few arcseconds around the 
non-thermal radio continuum source Sgr A*. They con- 
cluded that there is strong evidence for a ~2.5 10^ Mq 
central dark mass located within <0.015 pc of Sgr A* with 
a mass density of at least 6.5 10^ Mq pc~^. This high den- 
sity excludes the possibility that the central mass concen- 
tration is a central star cluster. Sgr A* is surrounded by 
a huge Hii region Sgr A West with a siz(^ of 2.1x2.9 pc, 
which was first observed by Ekers et al. (1975). The ion- 
ized gas in this region forms a spiral pattern (see e.g. Lo 
& Claussen 1983, Lacy et al. 1991) and is therefore called 
the Minispiral. The western side of this feature represents 
the inner edge of a large ring of molecular and atomic 
gas, the Circumnuclear Disk ( CND ) which extends from 
~2 to ~7 pc from the centre. It was observed by several 
authors: Gatley et al. (1986) (H2), Serabyn et al. (1986) 
(CO,CS), Giisten et al. (1987) (HCN), DePoy et al. (1989) 
(H2), Sutton et al. (1990) (CO), Jackson et al. (1993) (Oi, 
HCN), Marr et al. (1993) (HCN), Coil & Ho (1999, 2000) 
(NH3), and Wright et al. (2001) (HCN). The main results 
of their investigations are: 

— The CND has a hydrogen mass of a few lO'* M0, 

— the disk is very clumpy with an estimated volume fill- 
ing factor of <i>v ~ 0.01, 

— the clumps have masses of ~30 M©, sizes of ^^0.1 pc, 
and temperatures >100 K. 

It is now clear that the CND is interacting with the sur- 
rounding molecular clouds, but it is still a matter of debate 
where these connections are located (Zylka et al. 1990, 
Coil & Ho 1999 2000, Wright et al. 2001). 

Krolik & Begelman (1988) constructed a clumpy disk 
model for AGNs where cloud-cloud collisions arc respon- 
sible for the energy and momentum transport and thus for 
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the viscosity. Shlosman et al. (1990) discussed the possi- 
bility of fueling an AGN by the means of a cloudy disk. 
Sanders (1998) investigated the accretion process in the 
Galactic Centre on the basis of 'sticky particle' calcula- 
tions. He concluded that the gas features observed within 
10 pc from the Galactic Centre can be understood in terms 
of tidal capture and disruption of gas clouds on low angu- 
lar momentum orbits in a potential containing a point 
mass. The inner edge of the CND is explained by the 
formation of a 'dispersion ring', an asymmetric elliptical 
torus precessing counter to the direction of rotation. This 
feature can be maintained for many orbital periods. These 
calculations follow the orbital motions of 4000 particles, 
whereas the number of observed clouds in the CND is 
^--^500. Less clouds means less collisions and thus a lower 
viscosity (see e.g. Pringle 1981). This has influences on the 
way the dispersion ring forms. 

In a previous article, we constructed a self-consistent 
model to describe the physical and kinematical state of the 
CND (Vollmer & Duschl 2001a). The gas clouds were de- 
scribed as isothermal spheres embedded in an Hii region. 
The disk structure formed by the clouds was described 
by a quasi standard continuous accretion disk using ade- 
quately averaged parameters of the discrete cloud model. 
We succeeded in reproducing observed quantities as the 
cloud mass, the cloud radius, the electron density at the 
cloud edge, the electron density of the Hii region, the emis- 
sion measure, the disk height, and the dispersion velocity. 
A major result was that the coUisional time scale for one 
cloud is several Myr. The isolated CND might thus be 
much longer lived than previously assumed (Giisten et al. 
1987). 

In Vollmer & Duschl (2001b) we pointed out a possibil- 
ity to create the inner edge of the CND: as the clouds ap- 
proach the Galactic Centre only those with densities high 
enough to resist tidal disruption can survive. At a certain 
distance from the Galactic Centre, the clouds with high 
enough densities to survive become too heavy and will 
fragment and/or collapse. This critical distance depends 
on the radiation field because it determines the radius of 
the clouds. Assuming realistic conditions for the Galactic 
Centre, we ended up with an inner edge at ~2 pc. 
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The present article is devoted to the dynamics of the 
gas clouds in the CND. With the help of an improved 
numerical 3D coUisional cloud model we investigate the 
collisional time scales of the CND, which is assumed not 
to interact with the surrounding molecular clouds. A more 
realistic way to model the Galactic Centre region is to in- 
clude the molecular clouds, whose projected distances to 
the CND are small (see e.g. Mezger et al. 1989, Zylka et 
al. 1990). In the second part of this article we take this 
environment into account in our model and deduce more 
realistic collision rates for a gas clouds falling into the 
Galactic Centre with an already existing CND. A possi- 
ble scenario for the star formation history in the Galactic 
Centre is presented. 

2. The model 

The numerical simulations must take into account that the 
CND consists of several hundred gas clouds. The present 
code treats each cloud as a point with a given mass and 
radius. The mass-radius relation for the clouds might be 
given by the virial theorem 

rci = 0.478 X Mei[M0]/rci[K] pc (1) 
where is the cloud radius, Md the cloud mass, and 
Tci the cloud temperature. The temperature distribution 
with respect to the distance from the Galactic Centre is 
assumed to be 

Tci(i?) = 224//RH K (2) 
where R is the distance of the cloud from the Galactic 
Centre. This distribution is consistent with the observed 
UV radiation field (VoUmer & Duschl 2001a) and the de- 
duced gas temperature (see e.g. Glisten et al. 1987, Sutton 
et al. 1990). On the other hand, Vollmer & Duschl (2001a) 
have shown that the clouds in the CND have a constant 
radius ~ 0.05 pc if one takes the ionization front due 
to the external UV radiation field into account. In the 
present study we use both mass-radius relations. 

We follow the orbits of these clouds in the three di- 
mensional gravitational potential. The radial distribution 
of the total enclosed mass M{R) is given by 
M[R) = M + MoRi (3) 
where M = 3 10^ M0 is the mass of the central black hole, 
and Mo = 1.610^ Mo/pc* describes the mass distribu- 
tion of the stellar content. This is close to the findings of 
Eckart & Genzel (1996). When orbiting around the Galac- 
tic Centre, the clouds can have inelastic collisions. For the 
search of the next neighbour of a cloud a Barnes & Hut 
(1986) treecode is used. During these collisions clouds can 
exchange mass and larger clouds can grow through coales- 
cence. The result of a collision can be one (coalescence), 
two (mass exchange), or many fragments (fragmentation). 
In our model, we limit the fragmentation case to a maxi- 
mum of 3 fragments. Let the radius of the first cloud be ri , 
that of the second cloud r2. Let the impact parameter be 
6, the velocity of the fragment Uf, and the escape velocity 
Vesc- We follow the prescriptions of Wiegel (1994): 



— for ri — r2 < < ri + 7-2 : 
fragmentation 

— for 6 < ri — r2 and Vcsc > vi- 
mass exchange 

— for b < ri — r2 and Vcsc < "^f : 
coalescence 

In order to avoid the production of large number of too 
small clumps {M < 10 Mq), we treat a fragmentary col- 
lision, which would give rise to a cloud with M < 10 Mq, 
as a mass-exchange collision in adding the mass of the po- 
tential third cloud to that of the lightest colliding clouds. 
This procedure ensures that the arising cloud mass spec- 
trum has a maximum between 10 and 15 Mq. 

The integration of the ordinary differential equation is 
done with the Burlisch-Stoer method (Stoer & Burlisch 
1980) using a Richardson extrapolation and Stoermer's 
rule. This method advances a vector of dependent vari- 
ables y{x) from a point a; to a point a; -I- by a se- 
quence of n substeps. Thus, the initial timestep H is di- 
vided subsequently into n=2, 3, 4, etc. substeps. At the 
end the solution of y{x + H) is extrapolated and an er- 
ror can be estimated. The size of the timestep is adap- 
tive and linked to the estimated error of the extrapola- 
tion. This error is normalized by the values of the dis- 
tance covered during the last timestep and the velocity 
of each particle. The error level for acceptance of the 
extrapolated solution is a free parameter and has to be 
adapted to the physical problem treated. For a relative 
error level Ari/rnorm = Aui/unorm < 0.1, where i=l, 2, 3, 
''norm = 10~^ pc and Unorm = 5 10~^ kms~^ the Courant 
criterion is fulfilled for each timestep. 

The collisions are evaluated at each timestep h — H/n 
and only those, which appear for all sequences n, are taken 
into account. Therefore, the relative error level e is crucial 
for the value of the obtained collision rates. We adopt the 
strategy to chose the error level in a way to match the 
theoretical collision rates. 

3. The normalization of the collision rate 

With the fixed precision of the calculations it is possible 
to compare the model collision time for one cloud with the 
theoretical one: 

tcoU ^ (»T-clCrclUcl)~^ (4) 

where rid is the local cloud density, Cci = tt^i the cross 
section of the cloud, and Wci the cloud velocity. The com- 
parison of the mean collision time for one cloud are made 
with the simplest spatial configuration, namely a spheri- 
cal volume. The clouds are initially distributed uniformly 
in a sphere with a radius of 4 pc. The initial velocities are 
uniformly distributed between 20 kms~^ and 40 kms~^. 
The gravitational potential is due to the following mass 
distribution: M{R) = Mq Ri with Mq from Eq. (|). In 
order to check our numerics, we made one set of simula- 
tions where the potential collisions are counted, but not 
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executed. In this case there is only negligible angular mo- 
mentum transfer due to gravitational interactions between 
the clouds and the spatial distribution of the clouds stays 
approximately constant during the simulation. 

The influence of the mass-radius relation on the colli- 
sion rate is tested with two different relations: 



— RUNl: Tci (X McivR- In this case the cloud tempera- 
ture adjusts to the temperature given in Eq. (||). Since 
there is cloud heating due to cloud-cloud collisions, 
it is assumed that the dynamical time scale is much 
larger than the cooling time scale due to radiation. 

— RUN2: Tci = const.. In this case an external UV radi- 
ation field is assumed that creates an ionization front 
on the cloud surface which is directed to the Galactic 
Centre. The location of this ionization front does not 
depend on the clouds' distance to the Galactic Cen- 
tre (Vollmer & Duschl 2001a), i.e. the cloud radius is 
constant. 

The ratio of the model mean collision time of one cloud 
and the theoretical one for RUNl and RUN2 without ex- 
ecuted collisions is shown in Fig. 0. The theoretical mean 
collision time is given by 



i. 



thco 
coU 



{nc\(Jc\Vc\)'^ ~ ((nci)(a-ci)(uci)) ■ (5) 
The ratio for both simulations stays constant over 
100 Myr. It is -1 for RUNl and ~0.75 for RUN2. Thus, 
both simulations reproduce the theoretical collision time 
scale correctly within 25%. The difference between RUNl 
and RUN2 lies in the density gradient of the cloud dis- 
tribution. Clouds at the outskirts of the sphere are larger 
for RUNl. Therefore, the average theoretical collision time 
scale is smaller, and the ratio between the model and the 
theoretical collision time scales is larger. Since these clouds 
are located in a region with a low particle density and do 
not contribute much to the global model collision time 
scale, the ratio between the model and theoretical colli- 
sion time scale is larger. 

We then calculated another set of simulations with pro- 
ceeding cloud-cloud collisions. As the system evolves, the 
angular momentum transport through these collisions, i.e. 
the viscosity, makes the particle density increase in the 
centre. The initial and the final stage of such a simulation 
is shown in Fig. |[ The ratio between the model and theo- 
retical mean collision time for RUNl and RUN2 is shown 
in Fig. ^. The theoretical mean collision time is given by 
Eq. d). 
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Fig. 1. The ratio between the model and theoretical mean 
collision time without executed collisions. Left panel: 
RUNl. Right panel: RUN2. 
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Fig. 2. Right: initial state of a simulation with collisions. 
Left: final state of the simulation after 10^ yr. 



RUNl: During 6 Myr the ratio between the model and 
theoretical collision time scales is ~1. As the spatial 
cloud distribution changes to higher central densities, 
the number of collisions in the system is overestimated 
by Eq. (^, i.e. the coUisional time scale is underesti- 
mated. This makes the ratio between the model and 
the theoretical collision time scale increase. The critical 
cloud density for this increase is rid ~ 100 clouds/pc'^. 



— RUN2: We observe here the same trend as for RUNl. 
When the cloud distribution contracts, the number 
of collisions is theoretically underestimated. This ef- 
fect also starts at a critical cloud density of rid ^ 
100 clouds/pc^. 

We thus conclude that the mean collision time for one 
cloud, which is resulting from the model, reflects the real 
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Fig. 3. The ratio of the model mean colhsion time to 
the theoretical one with proceeding collisions. Left panel: 
RUNl. Right panel: RUN2. 



value or might underestimate it by ^--^25% if the cloud dis- 
tribution is not heavily concentrated. 

4. CND simulations as an isolated structure 

We have distributed 500 particles within a ring volume of 
constant height H=l pc and an inner radius of i?in=2 pc 
and an outer radius i?out=7 pc. The initial velocities were 
Keplerian using a gravitational potential made by the 
mass distribution described in Eq. (^). A turbulent ve- 
locity of 20% of the Keplerian value was added in an 
arbitrary direction. The system evolved then freely in 
the gravitational potential of a point mass (the central 
black hole) and an extended mass distribution (the cen- 
tral star cluster). VoUmer & Duschl (2001b) have shown 
that the clouds, which approach the Galactic Centre, be- 
come tidally disrupted if they are not dense enough to 
resist. Those clouds with a high enough central density to 
be stable at a radius of i?ciit ~2 pc will collapse and/or 
fragment. In any case, no cloud can survive at radii closer 
than i?crit- In our model we take this effect into account 
by removing clouds with R < -Rcrit- 
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Fig. 4. The evolution of the cloudy disk model without 
collisions. The elapsed time is shown at the top of each 
frame. The x-axis shows the distance of the clouds to the 
centre, the y-axis shows the vertical distance with respect 
to the disk plane. 
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Fig. 5. The evolution of the collision time scale for one 
cloud for the simulation shown in Fig. ^ 



As in Sect, ^we have made first a simulation counting 
the collisions without executing them. The cloud radius 
is Tci = const = 0.05 pc. Fig. |^ shows different states of 
evolution of such a disk-like structure. The timesteps are 
At = 16 Myr. The first plot shows the ring-like initial 
condition. After several Myr the disk reaches a state of 
equilibrium. The disk thickness at a distance of 5 pc is 
about 2 pc. The coUisional time scale for one cloud is 
shown in Fig. ^. It ranges between 2 and 3 Myr. 
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Fig. 6. The evolution of the cloudy disk model with colli- 
sions. The elapsed time is shown at the top of each frame. 
The X-axis shows the distance of the clouds to the centre, 
the y-axis shows the vertical distance with respect to the 
disk plane. 



In the next step we have made two simulations with 
collisions using the two different mass-radius relations 
RUNl and RUN2. Since the evolution of the spatial cloud 
distribution is quasi identical for RUNl and RUN2, we 
show only that of RUN2 in Fig. |. 

The observed thickness of the CND, i.e. 0.5 pc at a 
distance of 2 pc and 1.5 pc at a distance of 5 pc from the 
Galactic Centre (Glisten et al. 1987) lies between those of 
the model with and without collisions. The initial mass 
distribution 18 Mq < Afd < 38 Mg (Fig. g dotted line) 
smears out due to the inelastic collisions. The mass dis- 
tribution after 10 Myr is plotted as a solid line in Fig. ^. 
The maximum has shifted to smaller cloud masses and a 
high mass tail has been built. A part of these high mass 
clouds should have collapsed, but this mechanism is not 
included in our code. 

The resulting evolution of the model mean collision 
time can be seen for the mass-radius relation RUNl and 
RUN2 in Fig. ^. We observe an increasing mean collision 
time for the mass-radius relation RUNl (Fig. |8| left panel). 
This is, because the cloud distribution is shifted to smaller 
masses. Thus the clouds are smaller and the collision time 
scale, which is proportional to the square of the cloud 
radius, increases. For the case of a constant radius inde- 
pendent of the cloud mass (RUN2; Fig. || right panel) the 
collision rate increases, because the number of clouds and 
thus the cloud density increases. 
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Fig. 7. The mass spectrum of the clouds. Dotted line: 
initial mass spectrum. Solid line: mass spectrum after 
10 Myr. 



Since the cloud density is always smaller than fidit the 
results for the mean collision time scales is robust. The 
derived mean collision time scale for one cloud is of the 
order of 1-2 Myr. 

In VoUmer & Duschl (2001a) we gave a mean colli- 
sion time scale of ~10 Myr. The factor 5 between both 
approaches can be understood in the light of the analyti- 
cal model for turbulent, clumpy accretion disks (VoUmer 
& Beckert 2002). In this model, the viscosity is given 
as V — i?e^^uturb'driv, wlicre Wturb is the turbulent ve- 
locity dispersion and /driv is the driving wavelength for 
the turbulence. From the analytical model it follows that 
Zdiiv — H, where H is the disk height. Thus one obtains 
I' — Re~^vturhH . This viscosity prescription is also used 
in VoUmer & Duschl (2001a). The relation between the 
gas surface density E and the mass accretion rate M of 
the disk is (see, e.g., Pringle 1981) 

Stti^J: = M . (6) 
The collision time scale for one cloud is then tcoW = 
ReVl-^. In VoUmer & Duschl (2001a) we used Re = 1000, 
whereas a Reynolds number of Re — 20 has to be used 
in the framework of VoUmer & Beckert (2002), which is 
consistent with our dynamical results (Fig. ||). VoUmer & 
Duschl (2001a) assumed thus a relatively large Reynolds 
number. Consequently, they obtained a relatively large gas 
surface density (Eq. ||), i.e. a relatively large disk mass 
(Mgas ^ 1.510^ M0), about a factor of 10 larger than 
that what we have used here. 



5. Environmental effects 

In the previous Sections we have treated the CND as an 
isolated structure. In reality this is not the case. Observa- 
tion of molecular transitions (e.g. CO: Sutton et al. 1990; 
CO, CS: Serabyn et al. 1986, HON: CoU & Ho 1999, 2000, 
Wright et al. 2001) and mm observations (e.g. Mezger 
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Fig. 8. Model mean collision times for the mass-radius 
relation Left panel: RUNl. Right panel: RUN2. 



et al. 1989, Dent et al. 1993) have shown that there are 
several molecular cloud complexes in the vicinity of the 
Galactic Centre. Sgr A East Core, a compact giant molec- 
ular cloud with a gas mass of several 10^ M0 and the 
giant molecular clouds M-0. 13-0.08 and M-0.02-0.07 form 
the Sgr A cloud complex. We will compare our simula- 
tions to the observational results of Zylka et al. (1990), 
because their distinction of different cloud complexes is 
made on kinematical grounds. They obtained the follow- 
ing results: M-0. 13-0. 08 has radial velocities in the range of 
~5-25kms~^, a total mass of ~3 10^ Mq and lies in front 
of the Galactic Centre. M-0.02-0.07 consists of two differ- 
ent features, (i) Sgr East Core with a mass of ~2 10^ Mq 
and (ii) a curved streamer of total mass ^10^ Mq and ve- 
locities between 25 and 65 kms~^, which lies also in front 
of the Galactic Centre. Thus, there is neutral gas with 
a mass of several 10^ Mq, which is presumably located 
within the inner 50 pc around the Galactic Centre (see, 
however, von Linden et al. 1993). In addition, there are 
hints that they are kinematically connected to the CND. 
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Fig. 9. Initial conditions for the cloud infall into the 
Galactic Centre. The axis correspond to: left: right as- 
cension, declination; right: LOS (z-axis), declination. 

5.1. Initial conditions 

We let evolve a model of the CND as described above 
with 500 particles and a total mass of ~1.510^ Mq dur- 
ing 10 Myr. We then added a spherical mass distribution 
of 1000 particles, a total mass of 3 10^ Mg, and a size of 
10 pc at a distance of ^30 pc. The initial values of the 
velocity of these clouds were 30% of the Keplerian value 
with an additional velocity dispersion of 20% of this value. 
Due to the sub-Keplerian velocity, the cloud falls into the 
Galactic Centre colliding with the CND. This initial con- 
figuration can be seen in Fig. ^. We have made two differ- 
ent sets of simulations: 

(i) a prograde encounter, i.e. the orbital angular momen- 
tum of the cloud is close to that of the disk; 

(ii) a retrograde encounter, i.e. the orbital angular mo- 
mentum of the cloud is opposite to that of the disk. Tech- 
nically, this is realized in mirroring the z-axis for the disk 
clouds: = —2;', vl = —vl, where i is the number of the 
cloud. 

Furthermore, we varied the fraction of kinetic energy 
that is radiated away during a cloud-cloud collision. This 
has important influences on the collision and mass accre- 
tion rate of the system. 



5.2. The model evolution 
5.2.1. Retrograde encounter 

The evolution of a system with no loss of kinetic en- 
ergy per collision can be seen in Fig. |l^. The timestep 
is At=1.2 Myr. The cloud hits the CND for the first 
time a.t t ~ 0.3 Myr. The second collision takes place at 
t ~ 0.8 Myr. The period of the orbit is t^rh ^0.7 Myr. 
During its orbiting the cloud is stretched due to the tidal 
shear. The streamer becomes a filament of ~40 pc within 
1 Myr. As the cloud collides for the first time with the 
CND, the latter is heavily damaged and shows an az- 
imuthally asymmetric cloud distribution. Meanwhile, a 
counter rotating streamer is building up at the outer part 
of the CND. After --3 Myr the CND becomes less and 
less prominent. After ~2 Myr the streamer crosses itself 
for the first time and has a "brezel" form at t ^^3 Myr. 
Later on, it begins to form a triangular feature. At the 
end of the simulation {t ~10 Myr), we observe this tri- 
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Fig. 10. The evolution of the cloud infall into the Galactic 
Centre as the observer vifould see it from the Earth. The 
elapsed time is plotted on the top of each frame. 



angular feature containing a ring-like rotating structure 
extending from ^^5 to ^^10 pc and a counter rotating core 
(the former CND). The simulation with 10% energy loss 
per collision shows a more pronounced ring-like structure 
and a less prominent rotating core. 

This evolution can be also described by the number of 
collisions and by the mass accretion rate. We define this 
mass accretion rate by the accumulated mass of clouds per 
timestep At, which have distances to the Galactic Centre 
less than 2 pc. This evolution is shown in Fig. |ll] for four 
different energy loss rates per collision: 0%, 1%, 10%, and 
an energy loss that is proportional to the fraction of the 
clouds' mass, which participates in the collision (Krolik & 
Begelman 1988). The evolution of the mass accretion rate 
and the collision rate are similar for all energy loss rates. 
During the first collision at t r^OA Myr the total accreted 
mass is ~ 3000 Mq within 0.1 Myr. Nevertheless, this is 
not the major accretion event. This happens between 2 
and 3 Myr, when the tidally stretched cloud (streamer) 
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Fig. 11. 1st column: mass accretion rate; 2nd column: 
number of collisions per timestep At; 3rd column: ac- 
creted mass per timestep At and accumulate accreted 
mass. Rows: no energy loss, 1%, 10%, mass dependent 
loss of kinetic energy per collision. 



is crossing itself. The total accreted mass in this second 
event is ^lO** Mq within 2 Myr. At later stages of the sim- 
ulation the mass accretion rate decreases steadily having 
still some small peaks. The comparison with the collision 
rate shows that the first accretion event coincides with 
the the first collision of the cloud with the CND. The ma- 
jor event coincides with the second to fourth collisions and 
with the first crossing of the streamer with itself. The later 
and smaller accretion events coincide with later collisions 
between the streamer and the CND. The peak mass ac- 
cretion rate at the beginning is ~3 10^^ M© yr^^ falling 
to several 10~^ M© yr~^ after an elapsed time of 6 Myr. 
The peak of the major accretion event is ~10~^ M© yr~^ 
decreasing steadily until the end of the simulation. 

The evolution of the simulation with 1% energy loss 
per collision is almost identical to that without energy 
loss. With an energy loss rate of 10% per collision the 
mass is accreted approximately 2 times faster than for 
the simulation without energy loss. After the first peak 
the mass accretion rate stays at a value > 0.01 Mgyr"^. 
Then it decreases strongly from 2 to 5 Myr and rises again 
at 9 Myr. The collision rate shows the same behaviour. 
Clouds collide more frequently during 1 and 2 Myr as for 
the simulation without energy loss. At t = 6 Myr the colli- 
sion rate begins to rise from 410^^ yr^^ to 1.5 10^'^ yr^^ 
at 9 Myr. This happens, because the cloud distribution 
of the infalling streamer symmetrizes and forms a second 
circumnuclear disk with a mass of ~ 1.5 10*^ M©. 
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Fig. 12. The evolution of the cloud infall into the Galactic 
Centre as the observer vifould see it from the Earth. The 
elapsed time is plotted on the top of each frame. 

5.2.2. Prograde encounter 

The evolution of the system with no loss of kinetic en- 
ergy per collision can be seen in Fig. The timestep 
is Ai=1.2 Myr. In difference to the retrograde encounter, 
the CND is not damaged by the infalling cloud. During the 
evolution of the system the cloud-cloud collisions change 
the angular momentum of the CND and that of the in- 
falling cloud. With no energy loss per collision this align- 
ment is minimum, i.e. at the end of the simulations the 
majority of the CND clouds have still their initial angular 
momentum. We observe at the end of the simulations a 
triangular structure with an outer ring-like structure and 
the CND inside. Thus, one can still recognize two struc- 
ture with different angular momenta. For the simulations 
with an energy loss greater than 5% the situation changes. 
At the end of the simulation there is only one kinematical 
entity, a new CND with an angular momentum between 
that of the CND and that of the infalling cloud. The evo- 
lution of the system is shown in Fig. O for 6 different 




^^^^^ 



Fig. 13. 1st column: mass accretion rate; 2nd column: 
number of collisions per timestep At; 3rd column: ac- 
creted mass per timestep At and accumulate accreted 
mass. Rows: no energy loss, 1%, 5%, 10%, 20%, mass de- 
pendent loss of kinetic energy per collision. 

energy loss rates per colhsion e: 0%, 1%, 5%, 10%, 20%, 
and a mass dependent energy loss. Whereas in the case of 
a retrograde encounter the general behaviour of the dif- 
ferent quantities were the same for different energy loss 
rates, in the case of a prograde encounter these evolutions 
change dramatically for different energy loss rates. In the 
case of energy loss rates e <1% the mass accretion rate 
shows two peaks of M '--^ 2 10~^ MQyi~^. After a third 
small peak the mass accretion rate drops to zero. In the 
case of e=5% there are also two peaks at the beginning 
followed by an increase between t ~2 Myr and 4 Myr 
{M ~ 210"'^ M0yr~^) and decreases then slowly. For 
e=10% there is a gap in the mass accretion rate after the 
first collision between the infalling cloud and the CND. 
At later timesteps the mass accretion rate stays almost 
constant at a value of M ^ 2 10^^ M© yr~^. The gap dis- 
appears for e=20%. After a large first maximum between 
0.5 and 1.5 Myr (M 710"^ Mgyr-i), the mass ac- 
cretion rate has a minimum at t ^ 2.5 Myr, rises again 
to 710""^ M0yr~^ at t ~4 Myr and decreases again. In 
the case of a mass dependent energy loss the mass accre- 
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tion rate decreases steadily from its initial maximum of 
- 10-2 M0 yr-i to M ~ lO'^ Mg yr"! after 10 Myr. 




4-1 1 1 1 1 1 1 1 1 1 1 ^ 


tl 1 1 1 1 1 1 1 1 1 1- 


4-1 1 1 1 1 1 1 1 1 1 1 ^ 


\ f-3a.-301) \ 


\ (-35.-318) \ 


\ f-e3.-4Sl) \ 









<-so.-iat)) 



-100 100 -mo 100 -loo o loo -loo o loo -loo o loo 




600 400 200 -200 -400 -600 

R.A, offset (arcsec) 




-500 -400 -200 200 400 

LOS offset (orcsec) 



Fig. 14. The Sgr A Radio and GMC Complex. Greyscale: 

A 1.3 mm with the IRAM MRT (Mezger et al. 1989, Zylka ^^g- ^eft panel: combined snapshots of the high mass 

et al. 1990). Contour: the synchrotron sheU Sgr A East simulation for the comparison with the Sgr A cloud com- 

observed at A 6 cm with the VLA (Yusef-Zadeh & Morris, P^^x- Right Panel: the same mass distribution but with 

1987). Frame panels: i3C0(2-l) IRAM spectra (Mezger et ^^^^ hne of sight (LOS) as the j/-axis. Smaller distances 

al 1996) nearer to the observer. 



The evolution of the collision rate also shows big differ- 
ences for different energy loss rates e. Whereas the peak 
(*coii ~ 2 10"'^ yr""^) due to the first encounter between 
the infalling cloud and the CND is always there, the sub- 
sequent behaviours are markedly different. For e=0% the 
collision rate begins to rise at t=3 Myr until t=6 Myr 
and then stays constant it~^ii 5 10~^ yr~^)- For e=5% 
the collision rate begins to rise at t ~1.5 Myr and stays 
constant between 3 and 10 Myr (i~Qij ~ 2 10~^ yr~^)- For 
higher e the collision rate even decreases at t ~ 4 Myr. 

6. Comparison with the Sgr A cloud complex 

We have constructed a simple toy model in order to com- 
pare the observed Sgr A cloud complex (Zylka et al. 1990) 



to a snapshot of our simulation (for the observational in- 
terpretation see Sect. ||). The Sgr A Radio and GMC Com- 
plex observed at A 1.3 mm with the IRAM MRT (Mezger 
et al. 1989, Zylka et al. 1990), the synchrotron shell Sgr A 
East observed at A 6 cm with the VLA (Yusef-Zadeh & 
Morris, 1987), and 1^00(2-1) IRAM spectra (Mezger et 
al. 1996) are shown in Fig. |l^. Since there seems to be 
at least two different kinematical features, we have added 
two snapshots of the tidally stretched cloud at two differ- 
ent times (0.6 Myr and 1.6 Myr) together with the CND. 
This has been done using a simulation with 5 times more 
clouds in order to have a particle resolution. The cloud 
cross sections were adapted to ensure the same mean colli- 
sion time scale as for the simulation with 1500 clouds. We 
use a prograde simulation, because it fits best observed 
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line-of-sight locations of the GMCs in the vicinity of the 
CND. At these early stages of evolution, the results do 
only marginally depend on the energy loss rate per colli- 
sion e. We therefore chose e=0. 

This composite snapshot can be seen in the left panel 
of Fig. It should be stressed here that the aim of this 
comparison is not to reproduce each feature in detail but 
to make a generic picture of the Sgr A region. Therefore, 
it is not troublesome that the cloud in the east of the 
Galactic Centre is not as close as M-0.02-0.07. However, 
the CND and the streamer in south-north direction con- 
taining M-0.13-0.08 are nicely reproduced. The mass dis- 
tribution in direction of the line of sight (LOS) is shown 
in the lower panel of Fig. The streamer in south-north 
direction is located in front of the Galactic Centre as well 
as the eastern cloud. 
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Fig. 16. The radial velocity field of Fig. |5[ Dashed lines 
represent negative velocities, dotted lines positive veloc- 
ities. The contours are in the range of -120-120 kms~^ 
with a stepsize of 10 kms~^. 

The radial velocity field can be seen in Fig. 061 The 
velocities in the R.A.-Dec. plane are shown in Fig.|l^. The 
streamer in south-north direction is thus approaching the 
Galactic Centre from the front crossing the northern part 
of the CND. The eastern cloud is near the apocenter of 
its elliptical orbit, the north-western part having already 
changed the sign of the radial velocity falling onto the 
CND from the front. 

Zylka et al. (1998) show CS(3-2) spectral Hue data ob- 
served at the IRAM 30m Telescope covering the central 
21' X 20' with a spatial resolution of 17". Their galactic 
longitude-velocity and galactic latitude-velocity plots can 
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Fig. 17. The velocity field of Fig. [ij in the R.A.-Dec. 
plane. For clarity only each 10th velocity vector is plotted. 
The CND is rotating counter clockwise. Its outer rim is 
rotating clockwise. 



be directly compared with our simulations. We therefore 
made plots of exactly the same slices in AZ and A6 respec- 
tively. 

The galactic longitude-velocity plots of Zylka et al. 
(1998) are plotted in Fig. [l^. T he corresponding model 
diagrams can be seen in Fig. 19. The main features of the 



CS(3-2) {v-M) plots are reproduced by the model. Ob- 
viously, the CS(3-2) intensity does not match the column 
density of the model plots. This might be due to the fact 
that the CS(3-2) emission depends on the physical condi- 
tion out of which only one is the total column density. In 
Fig. |l^(a) the streamer in south-north direction is clearly 
visible. Nevertheless it appears more prominent in the cor- 
responding CS(3-2) plot (Fig.ga)). In Fig. |l|(b) and (c) 
the CND is prominent. We miss the observed prominent 
cloud complex at positive velocities and positive longi- 
tudes almost completely. The eastern model streamer can 
be seen in in the upper right part of Fig. |l^(a) and as an 
extension of the CND structure at positive velocities in 
Fig. |l^(b). Both features have a counterpart in Fig. |l8|(a) 
and (b). In contrast, the feature at negative velocities in 
Fig. p^(d), which is due to the eastern part of the east- 
ern streamer (Fig. 16), in Fig. |l^(d) has no counterpart 
in Fig. ^(d). This implies that the observations might 
correspond to a slightly later timestep in the evolution of 
the eastern model cloud and/or that its orbital parame- 
ters might be slightly different from that of the observed 
cloud. Nevertheless, our eastern model cloud represents a 
valuable first approach in modelling M-0.G2-G.07. 




Fig. 18. Representative galactic longitude-velocity plots {v-Al) of the CS(3-2) emission. The data have been averaged 
over selected ranges in galactic latitudes, Ab (Zylka et al. 1998). 
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Fig. 19. Corresponding galactic longitude-velocity model 
plots {v-td). 



The galactic latitude-velocity plots of Zylka et al. 
(1998) are plotted in Fig. pG. The corresponding model 
diagrams can be seen in Fig. |21|. 

We succeed in reproducing the feature at negative lat- 
itudes and positive velocities (Fig. |2l|(e) and (f)), which 
is caused by the southern streamer. The model velocities 



are offset by ^--^40 kms~^ from the observed ones in both 
diagrams. The eastern streamer, which can be seen at neg- 
ative latitudes in Fig. |l|(a) and (b) has only an observed 
counterpart in Fig. |2^(b) . This is due to the fact that the 
eastern model streamer is located more to the north than 
M-0.02-0.07. The additional CND features at negative lat- 
itudes and negative velocities (Fig. [2l](c) and (d)) might 
indicate that the position angle of the CND might differ 
by some degrees between the model and the real CND. 

7. Discussion 

In this work we investigated the dynamical behaviour of 
the CND in the Galactic Centre. In a first step we treated 
an isolated ring-like structure which contains 500 clouds 
with a cloud mass of ~30 Mq. Once in equilibrium, the 
disk shows the same height with respect to the radius as 
it is deduced from mm-observations (Giisten et al. 1987). 
We normalized the number of collisions per unit time with 
the help of a simple spherical cloud distribution. The re- 
sulting value for the CND simulation in equilibrium is 
^cou ~ 2 10^ yr. Therefore, the isolated CND can be longer 
lived as it was thought before. This collision time scale im- 
plies a mass accretion rate of several lO""* M0yr~^. We 
conclude that the current mass accretion rate at pc scale 
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Fig. 20. Representative galactic latitude-velocity plots {v-Ab) of the CS(3-2) emission. The data have been averaged 
over selected ranges in galactic longitude, Al (Zylka et al. 1998). 
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changes in the behaviour of these quantities. For small 
energy loss rates per collision (e <2%) only ~2 10^ M© 
are accreted within less than 2 Myr, whereas the collision 
rate; increases by a factor of ~5 between 2 and 5 Myr. 
For intermediate energy loss rates (5%< e <10%) we ob- 
serve an approximately constant mass accretion rate of 
M ~ 3 10"^ Mq yr^i with a gap {At -0.5 Myr) after the 
first encounter of the cloud with the CND. For e=5% there 
is second broader gap {At ~2 Myr) ~1 Myr after the first 
collision of the cloud with the CND. The collision rate be- 
tween the clouds stays constant (f^jj ~ 3 10~^ yr"^) over 
the whole evolution. For e >20% the initial value of the 
mass accretion rate is more than a factor of 2 higher than 
for e <20% and decreases steadily during the whole evolu- 
tion. The collision rate shows a prominent peak when the 
cloud collides for the first time with the CND, drops then 
to i^jj ^ 10~^ yr~^ and varies by a small amount during 
the further evolution. 



Fig. 21. Corresponding galactic latitude-velocity model 
plots {v-Ab). 

In a second step the environment of the CND was taken 
into account. We have constructed a dynamical model 
where a spherical cloud falls into the Galactic Centre 
onto a pre-existing CND. A retrograde and a progradc 
encounter were simulated with different kinetc energy loss 
rates e per collisions. In the case of a retrograde encounter 
the evolution of the mass accretion and cloud collision rate 
during the collision between the cloud and the CND are 
not very sensitive to e. We observe a first maximum during 
At —0.1 Myr followed a larger maximum {At >1 Myr). 
In the case of a prograde encounter we observe dramatic 



This strong dependence of the mass accretion and the 
collision rate on the energy loss rate per collision shows 
the importance of cooling mecanisms and the magnetic 
field, which both influence the energy loss rate. Whereas 
the Ol and Cll lines are the major contributors to radia- 
tive energy losses in the shock front during the collision, 
the magnetic field retains kinetic energy in the clouds. Fol- 
lowing Krolik & Begelman (1988) the loss rate of kinetic 
energy during a cloud-cloud collision is given by e = / ^, 
where / is the fraction of cloud's mass participating in the 
collision and ^ is the degree of inelasticity. Krolik & Begel- 
man (1988) estimate / = 0.2, which is in good agreement 
with our findings for the simulations with a mass depen- 
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dent energy loss rate. The degree of inelasticity is given 

by 

where Ai is the Mach number and /3 is the ratio between 
gas and magnetic pressure. A Mach number Ai ^ 20 and 
/3 ~ 1 gives ^ = 0.95. Only for (3 < 0.1 more than 20% of 
the energy available in the collision is retained by the mag- 
netic field. The strength of the magnetic field in the CND 
is about 1 mG (Yusef-Zadeh et al. 1996), the cloud den- 
sities are of the order 10^ cm""^, and the gas temperature 
within the clouds is --150 K (Vollmer & Duschl 2001a). 
This leads to (3 of the order of 1. Therefore, we conclude 
that magnetic fields might not play a preponderant role 
for the evolution of the CND as long as Eq. holds. 

Vollmer & Duschl (2001a) estimate the energy loss rate 
in the shock during the collisions to be e =10%. Thus, 
most realistic simulations are those with e=10-20% and a 
mass dependent energy loss rate (lower panels in Figs, pd] 
and 13). For e >10% the evolutions of the mass accretion 
and cloud collision rates of a retrograde and a prograde 
encounter are very similar: there is a prominent peak after 
the first collision between the cloud and the CND followed 
by steady decrease. For e =10% there are measurable dif- 
ferences between a retrograde and a prograde encounter 
(see Sect. [5^). 

In a further step we have constructed the observed 
features of the Sgr A cloud complex by adding two snap- 
shots of our simulations. We can reproduce the observed 
topography in a satisfying way in the R.A.-Dec. plane and 
in galactic longitude/latitude-velocity planes. This means 
that the orbit of the recently infalling model cloud is pro- 
grade with respect to the CND and the gas motion in the 
galactic disk. The observed warp of the CND (Glisten et 
al. 1987) can be natu rally e xplained by a former prograde 
encounter (see Sect. 5.2.2 ). Based on our model, there 
is the possibility that one of the streamers has already 
passed the Galactic Centre less than 0.5 Myr years ago. 
During this former interaction of a cloud and the CND a 
small starburst might have taken place at a distance of 
several pc from the Galactic Centre. Since the front end 
of the tidally stretched cloud and the front of the east- 
ern streamer are just crossing the CND today, we should 
be able to find signs of a small current starburst in the 
CND. This might be the case for the observed Oi cloud 
(the Tongue, Jackson et al. 1993). In this case, the incom- 
ing streamer hits the CND at ~3-5 pc. The clouds lose 
angular momentum due to collisions and are contracted. 
The so infalling perturbed clouds fragment, collapse, and 
form stars. These new ionizing O and B stars are still em- 
bedded in the part of the cloud from which they formed. 
Ott et al. (1999) found evidence for ongoing star forma- 
tion in the Central Cavity. Especially, they identified stars 
embedded in extended dust clouds within the Minispiral. 
The UV radiation of the young hot stellar population can 
penetrate efficiently into the fragmented cloud. Therefore, 



the gas in this region is predominantly in atomic form and 
can not be observed in HCN. 

The remaining question is: how and when did the Hci 
star cluster form? Genzel et al. (1996) found that the Hei 
stars are rotating in a direction opposite to that of Galac- 
tic rotation. If the star cluster was formed through an 
CND - molecular cloud interaction and if the stars' an- 
gular momentum reflects that of their parent cloud, this 
would imply that the colliding cloud was on a retrograde 
orbit. Since we have argued that the energy mass loss rate 
e >10% the two last two rows of Fig. ^ are relevant. If 
we assume that star formation is triggered by cloud-cloud 
collisions, the majority of stars is formed during the first 
2 Myr after the first collision between the cloud and the 
CND. Assuming a cloud mass of 30 Mq this leads to a 
total mass of colliding clouds of ~6 10^ Mq within 2 Myr. 
On the other hand, if the starburst produces stars with 
a given IMF, we can calculate the initial total gas mass 
needed to form the observed Hei stars. Assuming an ini- 
tial mass function $ = CM ~^ with cut-offs at 0.2 Mq 
and 50 M© and normalized by the observed Heii stars: 

1-50 Mq 

AfHcii ^ / $(m)dAf = 25 X M"'=" ~ 500 Mq (8) 

JlO Mq 

we deduce total stellar mass produced by the starburst of 
several 10^ Mq. This would imply that 30-50% of the 
colliding gas mass is transformed into stars. Since the 
cloud orbit is retrograde with respect to the rotation of 
the CND, the colliding clouds lose efficiently angular mo- 
mentum and thus fall into the Galactic Centre. Therefore, 
we expect that the stars are formed near the inner edge of 
the CND. They need ~7 Myr to build a centrally peaked 
spatial distribution (Vollmer & Duschl 2001b). After this 
time the initial cloud has the appearance of several stream- 
ers (see Fig. |l^ at t=7 Myr) and can not be recognized 
any more as one single cloud. Nevertheless, we can imag- 
ine that a part of Sgr A East Core could be identified as 
this remaining feature. In order to give an example for 
the dynamical behaviour of the system after 10 Myr we 
plot the velocity field of our simulation after this time 
(Fig. ^) . The mass distribution and the velocity field are 
complicated with one side approaching the observing and 
the other one receding. The position and inclination angle 
of the orbit of the infalling cloud might not be the one ob- 
served. Since it is difficult to separate really distinct kine- 
matical features in the observations, the discussed model 
offers only one plausible explanation for the complicated 
kinematical structure of Sgr A East Core. 

Our model still contains great simplifications concern- 
ing the physics of the cloud-cloud collisions (see Sect. ||). 
The effects of these collisions on the cloud population 
might change, for example, with a better knowledge of 
the detailed radiative magneto-hydrodynamics which is 
clearly needed. 

The shortcomings concerning the large scale dynamics 
of our model can in principle be overcome by modifying 
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Fig. 22. Velocity field after 10 Myr. This configuration 
can be compared to the velocity field of the molecular 
cloud Sgr A East Core. The dotted lines represent posi- 
tive radial velocities. The dashed lines represent negative 
radial velocities. 

the initial conditions for the infalling cloud. But, since the 
parameter space for these initial conditions is very large, 
it is not the aim of this work to reproduce each feature 
in great detail but to show general effects. Nevertheless, 
we can conclude that the Sgr A cloud complex consists of 
more than one dynamically independent feature. 

For the distinction of the dynamically different features 
in the Sgr A cloud complex we followed Zylka et al. (1990) 
who made this distinction on observational grounds. It 
is still possible that interferences of these features exist. 
Only an analysis of a complete data cube of this region 
together with the model calculation can shed more light 
on this question. 

8. Conclusions 

We have made numerical simulations using a collisional 

N-body code where each particle has a given mass and 
radius. These clouds can have inelastic that are described 
schematically. We have first normalized the collision time 
scale in our code using a spherical configuration. Simu- 
lations of an isolated disk structure corresponding to the 
Circumnuclear Disk (CND) in the Galactic Centre yields 
a mean collision time scale of one cloud of tcoii ^ 2 Myr 
and a mass accretion rate of 10~^ yr~^ < < 
10-3 M0yr-i. 

Since two giant molecular clouds (CMC) are located 
closer than 50 pc from the CND, we take in a second step 



this environment of the CND into account. We made sim- 
ulations for a retrograde and a prograde encounter using 
different loss rates of kinetic energy during the collisions. 
We found that the influence of the energy loss rate on the 
mass accretion and cloud collision rates is strongest for a 
prograde encounter. We estimate the energy loss rate per 
collision of the clouds in the Galactic Centre to be of the 
order of 10-20% of the kinetic energy. 

The direct comparison of a 1.2 mm maps and CS(3- 
2) position-velocity diagrams of a region containing the 
CND and the GMCs with a model snapshot shows that 
the GMCs are on a prograde orbit with respect to the 
rotation of the CND. A former prograde encoimtcr of a 
infalling cloud with the CND can naturally explain its 
observed warp. Since both GMCs might begin to collide 
with the CND we expect small current star formation near 
the inner edge of the CND. 

Within our scenario of an encounter between an in- 
falling molecular cloud and a CND, we point out the pos- 
sibility that the Hen star cluster has been formed by a 
retrograde encounter of a cloud with the CND ~7 Myr 
ago. The cloud that formed the Hell star cluster has been 
destroyed by tidal forces and can presently no longer be 
distinguished as one single kinematical entity. 
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